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Lepton flavour universality test at the CERN NA62 experiment 

Evgueni Goudzovski a on behalf of the NA62 collaboration 

a School of Physics and Astronomy, University of Birmingham, 
Edgbaston, Birmingham, B15 2TT, United Kingdom 

A precision test of lepton universality by measurement of the ratio Rk of K e + v to K + — > n + v decay rates 
was performed using a sample of 59963 K + — > e + v candidates with 8.8% background contamination collected by 
the CERN NA62 experiment. The result R K = (2.486 ± 0.013) x 10" 5 is in agreement with the Standard Model 
expectation. 



1. Introduction 

Decays of pseudoscalar mesons to light leptons 
are suppressed in the Standard Model (SM) by 
angular momentum conservation. In particular, 
the SM width of P ± -> l±v decays (with P = 
n,K,D( s) ,B) is 
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where Gf is the Fermi constant, Mp and Mi are 
meson and lepton masses, fp is the decay con- 
stant, and V qq , is the corresponding CKM matrix 
element. 

Within the two Higgs doublet models (2HDM 
of type II) , including the minimal supersymmetric 
one, the charged Higgs boson (H ± ) exchange in- 
duces a tree- level contribution to (semi)leptonic 
decays proportional to the Yukawa couplings of 
quarks and leptons [T]. In P ± — > i^v, it can 
compete with the W /± exchange due to the he- 
licity suppression of the latter. At tree level, the 
H exchange contribution is lepton flavour inde- 
pendent, and for P = tt, K, B leads to [2] 
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Here M# is the _ff ± mass, tan/3 is the ratio of 
the two Higgs vacuum expectation values, and 
e ks 10~ 2 is an effective coupling. 

A plausible choice of parameters Mh = 
500 GeV/c 2 , tan/3 = 40 leads to - 30% relative 



suppression of B + — > t + v decays, and ~ 0.3% 
suppression of K + ,D+ — » i + v decays with re- 
spect to their SM rates. However, searches for 
new physics in the decay rates are hindered by 
the uncertainties of their SM predictions. In par- 
ticular, interpretation of the measurements of the 
ratio T(K+ -> n+v)/T(K+ -> tt^^+v) in terms 
of constraints on (Mh, tan/3) phase space is cur- 
rently limited by lattice QCD uncertainties [3]. 

On the other hand, the ratio of kaon leptonic 
decay widths Rk = T(K e 2)/T(K^2), where the 
notation K e2 is adopted for K + — > i + v decays, is 
sensitive to loop-induced lepton flavour violating 
(LFV) effects via the H ± exchange [J] : 
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where the mixing parameter between the 
perpartners of the right-handed leptons 
can reach ~ 10~ 3 . This can enhance Rk 
by 0(1%) relative without contradicting any 
presently known experimental constraints, in- 
cluding upper bounds on the LFV decays r — > eX 
with X — 7], 7, /i/i. 

Unlike the individual decay widths, the ra- 
tio Rk = T(K e 2)/T(K fJ/ 2) is precisely predicted 
within the SM due to cancellation of hadronic un- 
certainties [S]: 
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where <5i? QED = (-3.79 ± 0.04)% is an elec- 
tromagnetic correction due to the internal 
bremsstrahhmg (IB) process. 

The sensitivity to LFV and the precision of 
the SM prediction make Rk an excellent probe 
of lepton universality. The current world aver- 
age (based on final results only) R^ A = (2.490 ± 
0.030) x 10~ 5 is dominated by a recent KLOE 
result [B]. A precise measurement of Rk based 
on a part (40%) of the data sample collected by 
the CERN NA62 experiment in 2007 is reported 
here. This is an update of an earlier result ob- 
tained with the same data sample [7]. 

2. Beam, detector and data taking 

The NA48/2 experimental setup [8] has been 
used for the NA62 2007-08 data taking. Exper- 
imental conditions have been optimized for the 
K e 2/K^2 measurement. The beam line is de- 
signed to deliver simultaneous unseparated K + 
and K~ beams derived from the SPS 400 GeV/c 
primary protons. However, the muon sweeping 
system was optimized for the positive beam in 
2007, and the sample used for the present analy- 
sis was collected with the K + beam only. Posi- 
tively charged particles within a narrow momen- 
tum band of (74.0 ± 1.6) GeV/c are selected by 
an achromatic system of four dipole magnets with 
zero total deflection, pass through a muon sweep- 
ing system, and enter a fiducial decay volume con- 
tained in a 114 m long cylindrical vacuum tank. 

With about 1.8 x 10 12 primary protons incident 
on the target per SPS pulse of about 4.8 s dura- 
tion, the secondary beam flux at the entrance to 
the decay volume is 2.5 x 10 7 particles per pulse, 
of which 5% are kaons (K + ). The fraction of 
beam kaons decaying in the vacuum tank at nom- 
inal momentum is 18%. The transverse size of the 
beam within the decay volume is Sx = 5y = 7 mm 
(rms), and its angular divergence is negligible. 

Among the subdetectors located downstream 
the decay volume, a magnetic spectrometer, a 
plastic scintillator hodoscope (HOD) and a liquid 
krypton electromagnetic calorimeter (LKr) are 
principal for the measurement. The spectrome- 
ter, used to detect charged products of kaon de- 
cays, is composed of four drift chambers (DCHs) 



and a dipole magnet. The HOD producing fast 
trigger signals consists of two planes of strip- 
shaped counters. The LKr, used for particle iden- 
tification and as a veto, is an almost homoge- 
neous ionization chamber, 27Ao deep, segmented 
transversally into 13,248 cells (2x2 cm 2 each), 
and with no longitudinal segmentation. A beam 
pipe traversing the centres of the detectors allows 
undecayed beam particles and muons from decays 
of beam pions to continue their path in vacuum. 

A minimum bias trigger configuration has been 
employed. The K e 2 trigger condition consists of 
coincidence of hits in the two HOD planes (the Qi 
signal), loose lower and upper limits on DCH hit 
multiplicity (the 1-track signal), and LKr energy 
deposit (i?LKr) of at least 10 GeV. The AT m2 trig- 
ger condition requires a coincidence of the Q\ and 
1-track signals downscaled by a factor D = 150. 

3. Measurement strategy, event selection 

The analysis strategy is based on counting the 
numbers of reconstructed K e 2 and candi- 
dates collected concurrently. Thus the analysis 
does not rely on the absolute beam flux measure- 
ment, and several systematic effects (e.g. due to 
charged track reconstruction and Q\ trigger ef- 
ficiencies, time-dependent effects) cancel at first 
order. 

Due to the significant acceptance and back- 
ground dependence, the measurement is per- 
formed independently in 10 bins of lepton mo- 
mentum covering a range from 13 to 65 GeV/c. 
The first bin spans 7 GeV/c, while the others are 
5 GeV/c wide. The selection conditions have been 
optimized separately in each momentum bin. The 
ratio Rk in each bin is computed as 

j_ N(K e2 ) - N B (K e2 ) A(K^) 
K D' N(K^)-N b (K^)' A(K e2 ) X 

f e X e(K e2 ) / L Kr ' 

where N(Ki2) are the numbers of selected Kg2 
candidates (£ — e,/z), Nb[K(2) are the numbers 
of background events, A{K^2) / A{K e2 ) is the ge- 
ometric acceptance correction, fi are the lepton 
identification efficiencies, t(Kg2) are the trigger 
efficiencies, /l& is the global efficiency of the LKr 
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readout, and D = 150 is the Ku2 trigger down- 
scaling factor. 

To evaluate the acceptance correction and the 
geometric parts of the acceptances for back- 
ground processes, a detailed Monte Carlo (MC) 
simulation including beam line optics, full de- 
tector geometry and material description, mag- 
netic fields, local inefficiencies of DCH wires, non- 
working LKr cells (0.8% of channels) and tem- 
porarily masked LKr cells is used. Particle identi- 
fication, trigger and readout efficiencies are mea- 
sured directly from data. 

Charged particle tracks are reconstructed from 
hits and drift times in the spectrometer. Track 
momenta are evaluated using a detailed magnetic 
field map. Clusters of energy deposition in the 
LKr arc found by looking at the maxima in the 
digitized pulses from individual cells in both space 
and time, and accumulating the energy within a 
radius of 11 cm. Shower energies are corrected for 
energy outside the cluster boundary, energy lost 
in non- working cells, and cluster energy sharing. 

Due to the topological similarity of K e 2 and 
K^2 decays, a large part of the selection is com- 
mon for the two modes: (1) exactly one recon- 
structed particle of positive electric charge geo- 
metrically consistent with originating from a kaon 
decay is required; (2) extrapolated track impact 
points in the DCHs, HOD and LKr must be 
within their geometrical acceptances; (3) track 
momentum must be in the range (13; 65) GeV/c, 
where the lower limit assures the efficiency of 
the i?LKr > 10 GeV trigger condition; (4) no 
LKr energy deposition clusters with energy E > 
Eveto = 2 GeV and in time with the track are al- 
lowed unless they are consistent with being pro- 
duced by the track via direct energy deposition 
or bremsstrahlung; (5) the reconstructed decay 
vertex longitudinal position must be within the 
nominal decay volume; (6) distance between the 
charged track and the nominal kaon beam axis 
must be below 3.5 cm. 

The following two principal criteria are used 
to distinguish K e % from decays. Kinematic 
identification is based on constraining the recon- 
structed squared missing mass in positron (muon) 
hypothesis -Mf < M^ iss (£) = (P K - P £ ) 2 < Mf, 
where Pk and Pg are the four-momenta of the 
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Figure 1. (a) M^ liss (e) vs lepton momentum for 
reconstructed K e 2 and K^a decays; (b) E/p spec- 
tra of positrons and muons. 

kaon (defined as the average one monitored with 
K + — > 7r + 7r + 7r~ decays) and the lepton (under 
the e + or /i + mass hypothesis). The limits Mf 
and Aff have been optimized taking into account 
the resolution and backgrounds, and vary among 
lepton momentum bins in the ranges (0.013; 
0.016) and (0.010; 0.014) (GeV/c 2 ) 2 , respectively. 
Lepton identification is based on the ratio E/p of 
track energy deposition in the LKr to its momen- 
tum measured by the spectrometer. Tracks with 
(E/p) min < E/p < 1.1, where (E/p) mhl = 0.95 
for p > 25 GeV/c and (E/p) min = 0.9 other- 
wise, are identified as positrons. Tracks with 
E/p < 0.85 are identified as muons. 

4. Backgrounds 

Kinematic separation of K e 2 from decays 
is achievable at low lepton momentum only (p < 
30 GeV/c), as shown in Fig. la. At high lepton 
momentum, the decay with the muon mis- 
identified as positron (E/p > 0.95, as shown in 
Fig. lb) due to 'catastrophic' bremsstrahlung in 
or in front of the LKr is the largest background 
source. In order to measure the mis-identification 
probability P Me , a muon sample free from the typ- 
ical ~ 10~ 4 positron contamination due to /x — > e 
decays has been collected: a 9.2Xq thick lead (Pb) 
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wall covering ~ 20% of the geometric acceptance 
was installed in front of the LKr during a pe- 
riod of data taking. The positron component in a 
sample of muon candidates consistent with orig- 
inating from K^2 decays, traversing the Pb wall, 
with p > 30 GeV/c and E/p > 0.95 is suppressed 
to a negligible level (~ 10 -8 ) by energy losses in 
the Pb wall. 

However, the muon passage through the Pb 
wall affects the measured Pj e b via two princi- 
pal effects: 1) ionization energy loss in Pb de- 
creases P^e and dominates at low momentum; 2) 
bremsstrahlung in Pb increases P^ and domi- 
nates at high momentum. To evaluate the cor- 
responding correction factor /pb = P^ e /P^, a 
dedicated Geant4 based [9] MC simulation of 
muon propagation downstream the spectrometer 
involving all electromagnetic processes, includ- 
ing muon bremsstrahlung [10] . has been devel- 
oped. The relative systematic uncertainties on 
P^ e and Pj e b obtained by simulation are esti- 
mated to be 10%, mainly due to the simulation of 
cluster geometry and calibration (measured and 
simulated Pj b are shown in Fig. 2). However, the 
error affecting their ratio is significantly smaller 
(<5/pb//pb = 2%) due to partial cancellation of 
uncertainties. 

The K^2 background contamination has been 
computed to be (6.10±0.22)% using the measured 
Pj b corrected by the simulated /pb, and correct- 
ing for the correlation between the reconstructed 
M^ iss (e) and E/p. The uncertainty comes from 
the limited size of the data sample used to mea- 
sure Pj b (0.16%), the uncertainty Sf Ph (0.12%), 
and model-dependence of the M^ iss (e) vs E/p 
correlation (0.08%). 

The K^2 decay also contributes to background 
via the /i — > e decay in flight. Energetic forward 
daughter positrons compatible with K e 2 kinemat- 
ics and topology are suppressed by muon polar- 
ization effects [TT]. Radiative corrections to the 
muon decay [12] lead to a further ~ 10% rel- 
ative background suppression. The background 
contamination has been estimated to be (0.27 ± 
0.04)%. 

The structure-dependent (SD) K + — > e + v^ 
process [T3], not suppressed by angular momen- 
tum conservation (more specifically, its SD + corn- 
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Figure 2. Muon mis-identification probability 
with the Pb wall installed P^ b for (P/p)mm = 
0.95 vs momentum: measurement (markers with 
error bars) ; simulation with its uncertainty (solid 
and dashed lines). 

poncnt corresponding to positive photon helic- 
ity), represents a significant background source. 
A recent measurement of the K + — > e + u^j (SD + ) 
differential decay rate [5J has been used to evalu- 
ate the background contamination to be (1.15 ± 
0.17)%. The dominant uncertainty comes from 
the uncertainty on the rate of K + — > e + vj (SD + ) 
decay, which has been increased by a factor of 3 
with respect to that reported in [BJ, as suggested 
by a stability check of Rk with respect to a vari- 
ation of the P vc to limit. 

The beam halo background is induced by halo 
muons undergoing /i — e decays in the vac- 
uum tank, or being mis-identified as positrons. 
It has been measured directly by reconstruct- 
ing the candidates from a K~ data sam- 
ple collected with the K + beam (but not its 
halo) blocked, and a special data sample collected 
with both beams blocked. The control sample 
is normalized to the data in the region —0.3 < 
M mi ss (M) < -0.1 (GeV/c 2 ) 2 populated predom- 
inantly by beam halo events. The 'cross-talk' 
probability to reconstruct a K^ 2 candidate due 
to a K~ decay with e + emission (K~ — > Ti a D l^v 1 
K~ — > it~tt d , K~ — > l~ve + e~ , where de- 
notes the 7T° Dalitz decay ir° — > r )e + e~) is at 
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Figure 3. Reconstructed squared missing mass 
Af^ iss (e) distribution of the K e2 candidates com- 
pared with the sum of normalised MC signal and 
background components. 

Table 1 

Summary of backgrounds in the K e2 sample. 



Source 


N B /N(K e2 ) 




(6.10 ±0.22)% 


(p -> e) 


(0.27 ±0.04)% 


K e2l (SD+) 


(1.15 ±0.17)% 


Beam halo 


(1.14 ±0.06)% 


K e3 


(0.06 ±0.01)% 


K 27r 


(0.06 ±0.01)% 



Total background (8.78 ± 0.29)% 

the level of ~ 10~ 4 and is taken into account. 
The halo background contamination has been es- 
timated to be (1.14 ± 0.06)%, where the uncer- 
tainty comes from the limited size of the con- 
trol sample and the uncertainty of its normal- 
ization. The beam halo is the only significant 
background source in the K^ 2 sample, measured 
to be (0.38 ± 0.01)% with the same technique. 

The number of K( 2 candidates is N(K e2 ) = 
59, 963 (about four times the statistics collected 
by KLOE) and N(K m2 ) = 1-803 x 10 7 . The 
M^ iss (e) distributions of data events and back- 
grounds are presented in Fig. 3; backgrounds in 
the K e2 sample integrated over lepton momentum 
are summarized in Table 1. 



5. Systematic uncertainties 

The ratio of geometric acceptances 
A(K ^2) / 1 A{K e2 ) in each lepton momentum bin 
has been evaluated with a MC simulation. The 
radiative K + — > e + v^ (IB) process is simulated 
following [13] with higher order corrections ac- 
cording to 14 15 . Lepton tracking inefficiency 
due to interactions with the spectrometer mate- 
rial is included into the acceptance correction, 
and its simulation has been validated with the 
data. The main sources of systematic uncer- 
tainty of the acceptance correction are the lim- 
ited knowledge of beam profile and divergence, 
accidental activity, and the simulation of soft 
radiative photons. A separate uncertainty has 
been assigned due to the limited precision of the 
DCH alignment. 

A sample of ~ 4 x 10 7 positrons selected kine- 
matically from K + — > ir a e + i; decays collected 
concurrently with the main Kg 2 data set is used 
to calibrate the energy response of each LKr cell, 
and to study f e with respect to local position and 
time stability (in the kinematically limited mo- 
mentum range p < 50 GeV/c). A sample of elec- 
trons and positrons from the 4x 10 6 Kl — > ir^e^v 
decays collected during a special short (15h) run 
with a broad momentum band K\ beam allows 
the determination of f e in the whole analysis mo- 
mentum range. The measurements of f e have 
been performed in bins of lepton momentum; sep- 
arate measurements have been performed for sev- 
eral identified groups of LKr cells with higher lo- 
cal inefficiencies. The inefficiency averaged over 
the K e2 sample is 1 — f e = (0.73 ± 0.05)%, where 
the uncertainty takes into account the statistical 
precision and the small differences between K + 
and K® L measurements. 

The efficiency of the Q\ trigger condition has 
been measured using K^ 2 events triggered with 
a special control LKr signal: integrated over the 
sample, it is (1.4 ± 0.1)%. Owing to its ge- 
ometric uniformity, and the similarity of the K e2 
and distributions over the HOD plane, it 

mostly cancels between the K e2 and K^ 2 sam- 
ples, and the residual systematic bias on Rk is 
negligible. The inefficiency of the 1-track trigger 
for Kg 2 modes is negligible. The trigger efficiency 
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Figure 4. Measurements of Rk in lepton momen- 
tum bins, and the averaged Rk (indicated by a 
band). The uncertainties in momentum bins are 
partially correlated, and include statistical and 
systematic contributions. 

Table 2 

Summary of the uncertainties on Rk- 



Source 


5R K x 10 s 


Statistical 


0.011 


K^2 background 


0.005 


K + ->■ e+i/7 (SD + ) background 


0.004 


Beam halo background 


0.001 


Acceptance correction 


0.002 


Spectrometer alignment 


0.001 


Positron identification 


0.001 


1-track trigger efficiency 


0.002 



correction e{K ^2) / t{K e 2) is determined by the ef- 
ficiency e(i?LKr) of the LKr energy deposit trigger 
signal E'LKr > 10 GeV, which has been measured 
to be 1 - e(£ L Kr) = (0.41 ± 0.05)% in the first 
lepton momentum bin of (13; 20) GeV/c, and to 
be negligible in the other momentum bins. The 
corresponding uncertainty on Rk is negligible. 

Energetic photons not reconstructed in the LKr 
may initiate showers by interacting in the DCH 
or beam pipe material, which causes the DCH 
hit multiplicities to exceed the limits allowed by 
the 1-track trigger condition. This suppresses 
the K + — ¥ e + vj (SD + ) background by about 
10% relative (varying over the positron momen- 
tum). Evaluation of the 1-track inefficiency for 
K + — > e + wf (SD + ) partially relies on simulation; 



its uncertainty has been propagated into Rk- 

The global LKr readout inefficiency has been 
measured using an independent readout system 
to be 1 - / L Kr = (0-20 ± 0.03)% and stable in 
time. 

6. Result and conclusions 

The independent measurements of Rk in the 
10 lepton momentum bins and the average over 
the bins are displayed in Fig. 4. Extensive stabil- 
ity checks in bins of kinematic variables, against 
variation of selection criteria and analysis proce- 
dures have been performed. The uncertainties of 
the combined result are summarized in Table 2. 
The result is 

R K = (2.486 ± 0.011 stat . ± 0.008 syst .) X 10~ 5 
= (2.486 ± 0.013) x 10~ 5 . 

This is the most precise measurement to date; it 
is consistent with the SM expectation. 
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